Area dependence of interlayer tunneling in strongly correlated 
bilayer 2D electron systems at u T = 1 
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The area and perimeter dependence of the Josephson-like interlayer tunneling signature of the 
coherent vt = 1 quantum Hall phase in bilayer two-dimensional electron systems is examined. 
Electrostatic top gates of various sizes and shapes are used to locally define distinct vt = 1 regions 
in the same sample. Near the phase boundary with the incoherent vt = 1 state at large layer 
separation, our results demonstrate that the tunneling conductance in the coherent phase is closely 
proportional to the total area of the tunneling region. This implies that tunneling at vt — 1 is a 
bulk phenomenon in this regime. 
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Bilayer two-dimensional electron systems (2DESs) in 
a large perpendicular magnetic field B support an un- 
usual collective phase when the separation between the 
layers is sufficiently small, the temperature is sufficiently 
low, and the total density tit = ri\ + n-z of electrons 
in the system equals the degeneracy eB/h of a sin- 
gle spin-resolved Landau level created by the magnetic 
fields. Interlayer and intralayer correlations are of com- 
parable importance in this Landau level filling factor 
vt = v\ + V2 = nr/ieB/h) = 1 phase. The sys- 
tem exhibits a number of striking physical properties, 
including Josephson-like interlayer tunneling 2 and van- 
ishing Hall and longitudinal resistances^^ when equal 
electrical currents are driven in opposition (countcrflow) 
in the two layers. The electronic correlations responsi- 
ble for these properties may be understood in a number 
of equivalent languages, including that of itinerant fer- 
romagnetism and that of excitonic Bose condensation 6 . 
Indeed, the remarkable transport properties observed in 
counterflow are very suggestive of excitonic superfluidity, 
although experiments have so far always detected small 
amounts of residual dissipation. 

Interlayer tunneling provides particularly dramatic ev- 
idence that bilayer 2DESs at vt = 1 are quite different 
at large and small interlayer separations. When the ef- 
fective layer separation (defined as the ratio d/l of the 
center-to-center quantum well spacing d to the magnetic 
length I = (H/eB) 1 / 2 ) is large, the interlayer tunneling 
conductance dl/dV is strongly suppressed around zero 
interlayer voltage^. This suppression reflects the strongly 
correlated nature of single layer 2DESs at high magnetic 
field. In essence, an electron which is rapidly injected into 
a 2DES at v < 1 is completely uncorrelated with the ex- 
isting iV-particle system and thus creates a highly excited 
iV+l-particle state which can only slowly relax to equilib- 
rium. Similar considerations attend the rapid extraction 
of an electron from a 2DES at high magnetic field. Since 
interlayer correlations are negligible at large d/l, the net 
tunneling conductance spectrum of the bilayer is a simple 
convolution of the suppression effects in the individual 
layers. In contrast, at sufficiently small d/l the inter- 



layer tunneling conductance dl/dV at vt = 1 displays a 
strong and extremely sharp peak centered at V = 0. The 
sharp peak in dl/dV reflects a near-discontinuity in the 
tunneling I — V curves at V = 0. This is reminiscent of 
the dc Josephson effect, but it remains unclear how close 
the analogy is. 

The tunneling features observed at vt — 1 and small 
d/l are widely believed to reflect the development of 
spontaneous interlayer quantum phase coherence among 
the electrons in the bilayer^iSiliii 2 .^^^^. Electrons 
are just as strongly correlated with their neighbors in 
their own layer as they are with electrons in the oppo- 
site layer. Within the exciton condensate language, every 
electron in the bilayer is bound to a correlation hole in 
the opposite layer and there is no energy cost associated 
with transferring an electron from one layer to the other. 
Crudely speaking, the Coulombic energy penalties asso- 
ciated with injection and extraction of electrons in the 
tunneling process are cancelled by the excitonic attrac- 
tion of the electron to the hole in the final state. 

The purpose of the present paper is to report on an 
investigation of how interlayer tunneling in the coherent 
vt = 1 phase is distributed across the sample area. Via 
electrostatic gating we are able to study regions of differ- 
ent size and shape in the same sample. For effective layer 
separations d/l not too much smaller than the critical 
value, approximately (d/l) c w 1.91 as T — > in our sam- 
ples, we find that the zero bias tunneling conductance, 
G(0) = dl/dV at V = 0, is essentially proportional to 
the sample area. This suggests that tunneling at vt = 1 
is a bulk phenomenon in our samples, at least close to the 
critical layer separation. In spite of the intuitive nature of 
this result, it is not obvious a priori. Indeed, the strong 
similarity between bilayer 2DESs at vt — 1 and Joseph- 
son junctions suggests that tunneling currents may be 
restricted to micron-scale regions close to the source and 
drain contacts and not exist in the bulk of the sample 9 . 
We discuss this and other theories of interlayer tunnel- 
ing following the description of our experiment and its 
results. 

The bilayer 2DES samples used in this experiment are 
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GaAs/AlGaAs double quantum wells grown by molecular 
beam epitaxy (MBE). Two 18 nm GaAs quantum wells 
are separated by a 10 nm Al .gGao.iAs barrier. Modula- 
tion doping with Si populates the ground subband of each 
well with a 2DES of nominal density n.1,2 ~ 5.5 x 10 10 
cm~ 2 and low temperature mobility /1 « 10 6 cm~ 2 /V s. 
The tunnel splitting Ag^g between the lowest symmetric 
and antisymmetric states in the double well potential is 
estimated 2 *^ to be less than 100 /xK; this is roughly six 
orders of magnitude smaller than the mean Coulomb en- 
ergy e 2 /e£ at vt = 1. Independent electrical contact to 
the individual layers is achieved via a selective depletion 
method^, thus allowing the measurement of interlayer 
tunneling via standard lock-in techniques. 

In order to minimize systematic errors in the tunnel- 
ing conductances due to sample-to-sample variations, we 
employ a technique which allows multiple, independently 
controllable tunneling regions to be established on a sin- 
gle sample. Two such samples, A and B, (both taken 
from the same parent MBE wafer) have been studied. 
In both cases a large electrostatic gate is deposited on 
the back side of the sample. This gate allows for global 
control of the electron density in the lower 2D electron 
gas in the active region of the sample. Smaller gates, of 
varying sizes, are deposited on the sample's top surface. 
These gates allow for local control of the top 2DES den- 
sity. By adjusting the gates appropriately it is possible 
to establish a density balanced vt = 1/2+1/2 = 1 state 
underneath any one of the top gates while the remainder 
of the sample is in an unbalanced configuration. 
Since tunneling is strongly enhanced at vt = 1 for d/£ < 
(d/£) c , but strongly suppressed^ where vt ^ 1, it is not 
difficult to separate out the conductance of the vt = 1 
region. In fact, as we show below, this separation can 
also be done at zero magnetic field. In sample A, the 
active region of the device is a 200 /im wide bar, fully 
underlaid by the large back gate. This bar is crossed by 
four rectangular top gates, with lengths 100, 50, 20, and 
10 /ini. The top gates thus define four different regions 
with four different areas that are proportional to the top 
gate length. In sample B the mesa bar has an inner sec- 
tion of width 100 /im and two outer sections of width 200 
/im. A 100 //m top gate crosses the inner section while 
two 50 /im top gates cross the outer sections. This sam- 
ple allows the comparison of regions with the same area 
but different mesa widths and gate lengths. Schematic 
diagrams of samples A and B are shown in Figs, la and 
3a, respectively. 

At zero magnetic field momentum and energy con- 
serving tunneling between parallel 2DESs can only oc- 
cur when the subband energy levels in the two quantum 
wells line up2i. If the two 2DESs have equal electron 
densities this occurs at zero interlayer voltage, V = 0. 
The width of the tunnel resonance is set by the lifetime 
of the electrons in the 2DESs; at low temperatures this 
is dominated by impurity scattering 2 ^. When the bilayer 
system is density imbalanced the tunnel resonance moves 
to a finite interlayer voltage proportional to the density 
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FIG. 1: (color online) (a) Schematic diagram of sample A. 
White rectangle represents 200 /im-wide mesa containing the 
bilayer 2DES; ohmic contacts (not shown) to the individual 
2DES layers are attached near both ends of the mesa. Grey 
rectangle under the mesa is the global back gate; hatched 
rectangles crossing above the mesa are the various top gates 
under which tunneling is established, (b) Zero voltage tun- 
neling conductance G(0) versus top gate bias taken at B — 
and with an applied back gate bias of —15.18 V. Peaks corre- 
spond, in order of decreasing height, to the 100, 50, 20, and 10 
fim top gates, (c) Background-subtracted peak heights from 
(b) versus top gate length. The line fitted to this data has an 
^-intercept of 3 ± 2 /im. 



difference. We exploit this fact to allow for clean sepa- 
ration of the tunneling conductance under the top gates 
on our samples from the "background" tunneling coming 
from other portions of the device. Figure lb illustrates 
this with data from sample A. For the data shown, the 
back gate bias has been set to -15.18 V; this lowers the 
density in the bottom 2DES and thereby imbalances the 
bilayer. With the top gates set to zero bias a measure- 
ment of the net tunneling conductance at zero interlayer 
voltage, i.e. G(0), yields a relatively small background 
value; the main resonance having been displaced to a fi- 
nite interlayer voltage by the imbalance. However, if the 
bias voltage on one of the top gates is swept, a resonance 
appears when the density of the top 2DES under that 
gate matches the density of the lower 2DES. This is what 
Fig. lb displays, there being one curve for each of the 
four top gates. Each gate produces a clear resonance in 
G(0) when biased to about -0.17 V 23 . These resonances 
sit atop a background conductance of about 50 nS which 
arises from off-resonant tunneling from those regions of 
the device not under the swept top gate. 

Figure lc shows the background-subtracted peak tun- 
neling conductance for each of the four top gate traces 
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shown in Fig. lb. The background subtraction is nec- 
essarily slightly different for each top gate. The tunnel- 
ing conductance observed with all top gates grounded 
gives the amount of tunneling that comes from the en- 
tire tunneling region when imbalanced. When a bias is 
applied to one of the top gates to create a balanced re- 
gion underneath, the net area of the device which re- 
mains imbalanced is reduced by the area of the biased 
top gate. We use this fact to calculate the correct back- 
ground for each top gate. The corrected peak heights 
vs. top gate length are plotted in Fig. lc. The data are 
well-fitted by a straight line with a small x-intercept of 
approximately 3 ± 2 /im. For the data shown in Fig. 1, 
rii =7i2 = 3.3 x 10 10 cm~ 2 at the resonance; but this 
same linear dependence on gate length is observed at all 
matched densities ranging from 2.9 to 3.4 x 10 10 cm~ 2 . 
We therefore conclude that tunneling is proportional to 
area at B = 0, as expected. The positive x-intercept 
in the peak height versus top gate length relation may 
indicate that boundary effects are reducing the effective 
tunneling area. For example, fringe fields along the edge 
of the gated regions could produce a density gradient 
over length scales comparable to the distance between 
the upper 2DES and the top gate (0.5 /jm), thus reduc- 
ing the area in the upper 2DES that has electron density 
matching that of the lower layer. 

At fx = 1 the situation is both simpler and more 
complex than at B — 0. The strong suppression^ of 
tunneling in the unbalanced regions of the sample where 
vt 7^ 1 and the strong enhancement at small d/l of the 
tunneling in the balanced vt = 1 regions under the top 
gates are both helpful in keeping the background tun- 
neling conductance (typically ~1 nS) a small fraction of 
the peak conductance at Vt = 1 (typically ~10 - 1500 
nS). In contrast, it is much more difficult to ensure that 
the bilayer is density balanced at vt = 1 than it is at 
zero magnetic field. At vt = 1 the tunnel resonance re- 
mains centered at zero interlayer voltage, regardless of 
the individual layer densities, provided v\ -\-Vi = vt = 1- 
Furthermore, the robustness of the vt = 1 state against 
small antisymmetric shifts in the individual filling factors 
{v\ — > v\ + Sv, V2 — > V2 — 8v) complicates the determi- 
nation that a balanced vt = 1 system is present. In 
order to establish the gate voltages required for balanced 
vt — 1 state underneath any given top gate, we first use 
an empirical calibration of the gates^i to choose approx- 
imate top and back gate voltages (Vtop and Vback) which 
will create v\ « V2 ~ 1/2 in each layer at some chosen 
magnetic field. Since the empirical calibrations are rea- 
sonably good, this results in a clear tunneling peak at 
V = 0, thus proving that the system is at least close 
to the desired balanced i*r = 1 configuration. Next, 
the zero bias tunneling conductance G(0) is recorded as 
both the Vtop and Vback are scanned across a window in 
voltage space. Examination of the resulting data clearly 
reveals contours of constant total filling factor. Small 
refinements of the magnetic field and the gate voltages 
are then made to first maximize G(0) vs. B and then 
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FIG. 2: (color online) Tunneling vs. gate length at vt = 1 in 
sample A at d/l = 1.81. Gate layout as shown in Fig. 1. (a) 
Tunneling conductance resonances (dl/dV vs. V) at T = 85 
mK. Gate lengths 100, 50, 20, and 10 /jm; tallest to shortest, 
(b) Peak heights plotted vs. top gate length; open circles 60 
mK, solid circles 85 mK, squares 125 mK. Straight lines are 
linear least squares fits. 



to perfect the symmetry 2 ^, of the tunnel resonance shape 
(dl/dV vs. V). 

Figure 2a shows typical tunneling conductance reso- 
nances (dl/dV vs. V) at vt = 1 in sample A. For these 
data the temperature is T — 85 mK and the effective 
layer spacing is d/i = 1.81. The four traces (tallest to 
shortest) correspond to tunneling under the four differ- 
ent top gates, with lengths 100, 50, 20, and 10 /im. As 
explained above, there is negligible background tunnel- 
ing coming from other parts of the sample. Figure 2b 
shows the height of the tunneling peak at zero interlayer 
voltage, i.e. G(0), vs. top gate length at three temper- 
atures, 60, 85, and 125 mK, for this same effective layer 
spacing. The data reveal a clear linear dependence on 
gate length. A small ^-intercept, 5 ± 2 /jm is obtained 
from the linear least squares fits. Data very similar to 
that shown in Fig. 2 has been obtained at effective layer 
spacings ranging from d/i — 1.70 to 1.88 and at tempera- 
tures from T = 60 mK to 300 mK. In all cases the height 
of the tunneling peak scales linearly with gate length. 

At smaller d/l deviations from the linear dependence 
on gate length begin to appear. However, these devia- 
tions are most likely artifacts arising from the reduced 
sheet conductivity of the 2D systems at these lower den- 
sities and the concomitant enhanced tunneling conduc- 
tance. When the sheet conductivity of the 2DESs no 
longer greatly exceeds the tunneling conductance, volt- 
age drops develop within the 2D layers. The resonance 
in dl/dV becomes broadened and suppressed in height. 
Since this effect is more pronounced in larger tunneling 
areas than in smaller ones, the tunneling peak height 
becomes sub-linear in gate length. For this reason we 
restrict our attention here to effective layer spacings and 
temperatures relatively close to the phase boundary 26 
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FIG. 3: (color online) (a) Schematic diagram of the mesa and 
gates in the active region of sample B. (b) Zero voltage tun- 
neling conductance G(0) versus top gate bias taken at B — 
and with an applied back gate bias of —17.08 V. (c) Tunneling 
conductance resonances (dl /dV vs. V) at vt = 1 with d/t = 
1.64 and T = 61 mK. In (b) and (c) the solid trace (black) 
corresponds to the 100 /im square tunneling region and dot- 
ted trace (red) to the 200 x 50 /im rectangular region. The 
red trace in (c) has been multiplied by 1.086, as described in 
the text. 



separating the interlayer coherent and incoherent states 
at Vt = 1- In this regime we are confident that sheet 
resistance effects are negligible. 

The data in Fig. 2 strongly suggest that the interlayer 
tunneling conductance in the coherent vt = 1 state is 
proportional to the area of the tunneling region. How- 
ever, since in sample A it is the length of the gated regions 
that is varied while the width (200 /im) is constant, it re- 
mains possible that the tunneling conductance is simply 
proportional to gate length, not area. For this reason we 
now turn to the results obtained from sample B. This de- 
vice, illustrated schematically in Fig. 3a, contains three 
top-gated regions; two are 200 /im wide by 50 /im long 
rectangles and the third is a 100x100 /im square. These 
three regions all have the same area but differ by a factor 
of two in gate length and width. Since the two rectan- 
gular gates are identical in size and shape, we focus here 
on the comparison of the inner square top gate and one 
of these outer rectangular gates. Direct comparisons of 
tunneling in the two rectangles justify this simplification. 

Figure 3b shows typical tunneling conductance reso- 
nances in sample B at zero magnetic field. As with sam- 
ple A, at B = it is most informative to plot G(0), the 
tunneling conductance at zero interlayer voltage, versus 
top gate bias at fixed back gate bias, here Vb ac k — —17.08 
V. Although sample B was taken from the same MBE 



wafer as sample A, its tunneling conductance per unit 
area is about a factor of ~ 2 smaller. We attribute this 
difference to known wafer-scale variations in the thick- 
nesses of the various epilayers in the sample^ 7 -. More 
significantly, Fig. 3b shows that the peak tunneling con- 
ductance of the 100 /im square region is about 10 per- 
cent larger, and the resonance width about 10 percent 
smaller, than that in the 200 x 50 /im rectangle. While 
these differences might be due to differing degrees of den- 
sity inhomogeneity in the two regions, their exact origin 
is unknown. In the following, when comparing the tun- 
neling conductances at Vt = 1 in the square and rectan- 
gular regions of sample B, we correct for this small effect 
by multiplying the tunneling data from the rectangular 
region by the appropriate factor deduced from the zero 
field data. Figure 3c illustrates the results of this com- 
parison of tunneling conductance spectra (dl /dV vs. V) 
at Vt = 1 for d/t = 1.64 and T — 60 mK; there is es- 
sentially no difference in the spectra from these two dif- 
ferently shaped regions. Very similar results have been 
found at densities corresponding to 1.60 < d/i < 1.79 
and temperatures from T = 60 to 300 mK. 

Taken together, the data presented in Figs. 2 and 3 
provide strong evidence that the enhanced tunneling con- 
ductance characteristic of the coherent vt = 1 bilayer 
state is proportional to the area of the tunneling region, 
at least over a limited range of d/t and temperature near 
the phase boundary. The possibility that the conduc- 
tance is instead proportional to the perimeter of the tun- 
neling region, or simply the length of that region along 
the mesa boundary is not consistent with the present 
findings. 

The intuitive character of our findings belies the subtle 
nature of tunneling in the coherent bilayer vt = 1 state. 
The long- wavelength effective Hamiltonian for this many- 
electron state is closely analogous to that of a Josephson 
junction. As pointed out by Stern, et al£ this anal- 
ogy leads to a sine-Gordon equation for the collective 
phase variable <p representing the orientation of the pseu- 
dospin magnetization in the ferromagnetically ordered 
stat o 17 ' 18 . In an ideal, disorder-free sample this equa- 
tion leads to the conclusion that interlayer tunnel cur- 
rents (proportional to sin0) are restricted to narrow re- 
gions near where the current enters and leaves the bilayer 
system^ 8 .. The width Aj of these regions is determined by 
the pseudospin stiffness and the tunnel splitting As^s, 
and is estimated to be in the few /im range. In contrast to 
the results presented here, this scenario would not result 
in a tunneling conductance proportional to area3^. 

More realistic models of the distribution of tunneling in 
the coherent vt = 1 state have recently been advanced by 
Fertig and Murthy 30 . and by Rossi, et al^. Both models 
attempt to incorporate the crucial effects of disorder on 
the tunneling. In Fertig and Murthy's picture 3 ^ fluctua- 
tions in the 2DES density lead to a complex network of 
channels and nodes, covering the entire sample, in which 
the coherent vt = 1 state exists. Tunneling occurs at 
the nodes of the network. Since the disorder length scale 



5 



is expected to be small (< 1 /im, set by the separation 
between the dopant layers and the 2DESs), the number 
of nodes, and therefore the net tunneling conductance, 
is predicted to be proportional to sample area in large 
devices, in agreement with our results. In the model of 
Rossi, et alM- disorder along the physical edge of the 
sample is expected to dominate when the coherent phase 
is well-developed (i.e. at small d/t and very low temper- 
ature). In this regime they predict that the tunneling 
conductance will be proportional to the length along the 
sample edge, not the total area. Although this is in con- 
flict with our observations, we stress that our results are 
confined to relatively large d/£, in fairly close proximity 



to the phase boundary. 

In summary, we have found strong evidence that inter- 
layer tunneling in the coherent bilayer vt = 1 phase near 
the critical layer separation is proportional to the area of 
the tunneling region, with some evidence for small edge 
corrections. We therefore conclude that tunneling in the 
coherent vt = 1 state, much like tunneling at zero mag- 
netic field, is a bulk phenomenon in our samples. 
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